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Executive Summary 


The National Weather Service Forecast Office in Melbourne, FL (NWS MLB) is responsible for providing 
meteorological support to state and county emergency management agencies across East Central Florida in the event 
of incidents involving the significant release of harmful chemicals, radiation, and smoke from fires. NWS MLB 
currently relies on a PC-based version of the National Ocean and Atmospheric Administration Hybrid Single- 
Particle Lagrangian Integrated Trajectory (HYSPLIT) model to provide such guidance. HYSPLIT provides a 
complete system for computing atmospheric trajectories, complex dispersion, and concentration simulations using 
either puff or Lagrangian particles. The NWS Spaceflight Meteorology Group (SMG) is responsible for all weather 
forecasts for Space Shuttle landings including cumulus cloud development and visibility issues from smoke and 
high-level trajectory forecasts for thunderstorm anvils. SMG has expressed interest in using HYSPLIT to provide 
trajectory forecasts in support of Space Shuttle landings. 

The PC-based version of HYSPLIT is run through a graphical user interface (GUI) for ease in setting up a 
trajectory, concentration, or deposition guidance. While the PC-based GUI is convenient for making single-use 
trajectory and concentration simulations it is difficult for forecasters to run and manage in an operational 
environment. To alleviate the difficulty in providing scheduled real-time trajectory and concentration guidance, the 
Applied Meteorology Unit (AMU) was tasked to configure a Linux version of HYSPLIT that ingests input from 
several National Center for Environmental Prediction (NCEP) model products and the Weather Research and 
Forecasting (WRF) model. These objectives were met successfully. 

The AMU configured HYSPLIT to provide forecasters with automated guidance from a number of source 
locations using meteorological input from the NCEP North American Mesoscale (NAM) and the Rapid Update 
Cycle (RUC) models. For the NAM guidance, FIYSPLIT was configured to run four times a day, from a variety of 
user defined source locations and forecast intervals. For shorter-range guidance, the AMU configured HYSPLIT to 
run eight times per day using the operational RUC products from various source locations and forecast times. The 
scheduled guidance is controlled through cron jobs that automatically download the NCEP model products, convert 
the meteorological grids into HYSPLIT binary format, run the model, and post-process the data to create output 
graphics, including an option to display the output in Google Earth™ or within Geographic Information System 
(GIS) applications. 

NWS MLB and SMG run locally-configured versions of the WRF model within the Environmental Modeling 
System (EMS). The WRF EMS model configurations are initialized with local observations and are run at finer 
spatial and vertical resolution than provided by the NCEP model guidance. WRF can provide input into HYSPLIT 
that may better predict the unique mesoscale environment in Florida and southern California than with the NCEP 
guidance. Therefore, the AMU was tasked to configure a utility program to convert EMS output into HYSPLIT 
binary format. Since SMG uses the Advanced Research WRF (ARW) within the EMS package, a converter recently 
made available in HYSPLIT could be used without modifications. On the other hand, NWS MLB is running the 
Nonhydrostatic Mesoscale Model (NMM) version of WRF within the EMS and no such conversion program for 
HYSPLIT exists. The AMU configured a conversion program for the NMM output by modifying existing source 
code used to convert the NAM model output into HYSPLIT format. The end product allows forecasters to 
automatically provide trajectory and concentration forecasts on a scheduled basis using either the NCEP products or 
locally configured versions of WRF and, therefore, provide timely information on hazardous conditions to their 
customers. 
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1. Introduction 


The National Weather Service Forecast Office in Melbourne, FL (NWS MLB) is responsible for providing 
meteorological support to state and county emergency management agencies across East-Central Florida in the event 
of incidents involving the significant release of harmful chemicals, radiation, and smoke from fires and/or toxic 
plumes into the atmosphere. Accurate and timely guidance is critical for decision makers charged with protecting the 
health and well-being of populations at risk. Information that can describe the geographic extent of areas possibly 
affected by a hazardous release, as well as to indicate locations of primary concern, offers better opportunity for 
prompt and decisive action. NWS MLB uses the National Oceanic and Atmospheric Administration (NOAA) 
Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model to provide trajectory, concentration, and 
deposition guidance during such events. 

Forecasters at the NWS Spaceflight Meteorology Group (SMG) in Houston, TX have also expressed interest in 
using HYSPLIT to assist with weather Flight Rules (FR) during Space Shuttle landings at the Kennedy Space Center 
in Florida and Edwards Air Force Base in California. In particular, SMG would provide low and mid-level 
HYSPLIT trajectory forecasts for cumulus clouds associated with smoke plumes, and high-level trajectory forecasts 
for thunderstorm anvils. Another potential benefit for both NWS MLB and SMG is using HYSPLIT concentration 
and deposition forecasts in dense fog situations related to smoke plumes from fires. 

For expedience in delivering dispersion forecasts in the diversity of operational situations, NWS MLB had 
become increasingly reliant on the PC-based version of HYSPLIT accessed through a graphical user interface 
(GUI). While the PC-based GUI offers unique advantages when compared to running the model from a command 
line, it is still difficult for forecasters to run and manage in an operational environment. Therefore, the first goal of 
this task for the AMU was to obtain and configure a Linux version of HYSPLIT that ingests the National Centers for 
Environmental Prediction (NCEP) guidance, such as the North American Mesoscale (NAM) and the Rapid Update 
Cycle (RUC) models. The AMU configured scripts to automatically download the NCEP model products, convert 
the meteorological grids into HYSPLIT binary format, run the model from several pre-selected latitude/longitude 
sites, and post-process the data to create output graphics. 

NWS MLB and SMG run locally configured versions of the Weather Research and Forecasting (WRF) model 
contained within the Environmental Modeling System (EMS; Rozumalski 2006). These locally configured versions 
of WRF are run at finer spatial and vertical resolution than provided by the NCEP model guidance. The WRF 
configurations are initialized with their own local data assimilation systems and can provide input into HYSPLIT 
that may better predict the unique mesoscale environment in Florida and southern California than the coarser 
resolution NCEP guidance. Therefore, the second goal of this task was to develop a utility to convert local WRF 
output into HYSPLIT binary format. The AMU created a conversion program for the WRF output by modifying 
existing source code used to convert the NAM products. This allows forecasters to automatically provide trajectory 
and concentration guidance on a scheduled basis using either NCEP products or a locally configured WRF model 
and, therefore, provide timely information on hazardous conditions to their customers. 

This report is organized as follows. Section 2 provides a brief technical description of the HYSPLIT model. 
Section 3 gives an overview of the local HYSPLIT model configuration that uses the operational NCEP model 
guidance for NWS MLB. Section 4 provides a description of the local HYSPLIT configuration using the WRF 
model contained within the EMS run at both NWS MLB and SMG. A summary of the project is given in Section 5. 
Finally, a user guide that describes the model structure and description of the various HYSPLIT variables is given in 
Appendix A. 
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2. HYSPLIT Model Description 

The HYSPLIT model has been developed over many years at the NO A A Air Resources Laboratory (ARL). The 
model is a complete system for computing atmospheric trajectories, complex dispersion and simulations using either 
puff or Lagrangian particles (Draxler 2007). HYSPLIT applications include tracking and forecasting the release of 
pollutant material, air-parcel trajectories, volcanic ash, and smoke from wildfires. The model is used operationally 
by NOAA through various forecasting agencies such as NCEP and local NWS offices. HSYPLIT uses gridded 
meteorological data to drive the trajectories and concentrations either forwards or backwards in time at regular time 
intervals. The software is written mainly in Fortran with several modules in C code. A complete description of the 
model, code, and software structure can be found in Draxler (1997, 2007). 

Dispersion of a pollutant is calculated by assuming either puffs or particles are released into the wind field. In 
the puff model, puffs expand until they exceed the size of a pre-defined meteorological grid cell and then split into 
several new puffs, each with a share of the pollutant mass. In the particle model, a fixed number of particles are 
moved by the mean horizontal and vertical wind field including a turbulent component. For default simulations, 
HYSPLIT assumes a puff distribution in the horizontal and particle dispersion in the vertical direction. This allows 
for greater accuracy of the vertical dispersion parameterization combined with the advantage of having an ever- 
expanding number of particles to represent the pollutant distribution. HYSPLIT trajectories are calculated by 
running the model without dispersion, thus computing the advection of a single pollutant particle within the mean 
wind. For both trajectory and concentration computations, the model can be run forward or backward in time in 
order to locate an unidentified pollutant source location. 

The gridded meteorological data fields required to run HYSPLIT can be obtained from archived data sources or 
forecast model output available at NCEP. Fortran preprocessors are required to convert the model output fields, 
usually in the Gridded Binary (GRIB) format, into HYSPLIT’ s own binary format, hereafter referred to as ARL 
format. GRIB is a mathematically concise data format commonly used to store historical and forecast weather data. 
The ARL data structure is a compressed binary direct-access format where multiple time-periods can be contained in 
each file for quick read access at model run time. Each time period contains a separate index record that includes 
grid definition, quality control (QC) flags, variable identification and vertical level information. The model can 
accept gridded meteorological data from multiple grid configurations including polar, global latitude/longitude such 
as in the Global Forecast System (GFS) model, and Lambert conformal such as in the NAM model. During 
conversion of the GRIB data records, the input meteorological fields are interpolated to an internal sub-grid chosen 
by the user to reduce memory requirements and increase computational speed. In addition, HYSPLIT calculations 
can be performed on nested multiple grids, usually specified from fine to coarse resolution, which allows users the 
flexibility of using their own locally configured model. 

The HYSPLIT model can be run interactively on the internet through the NOAA website at 
http :// www.arl .noaa. go v/HYSPLIT.php with appropriate permissions, or the source code can be downloaded and 
built on a Linux system. The website version has been configured with some limitations to avoid computational 
saturation of the NOAA server, but the PC or Linux version of the code has no such restrictions. The model includes 
a Tool Command Language (Tcl/Tk) GUI for ease in setting up a trajectory, concentration, or deposition simulation. 
The main GUI includes capabilities to download model input data from a defined File Transfer Protocol (FTP) 
server, convert the data to ARL format, run the trajectory or concentration model, and post-process the output 
graphically. 

The HYSPLIT post-processing programs include several graphical packages used to generate multiple color or 
black and white postscript images. Recent additions to the software have allowed users to output graphics directly to 
Google Earth™ format. This option creates a compressed Google Earth .kml file that allows for animated four- 
dimensional images. HYSPLIT includes the capability of displaying concentration or trajectory guidance on detailed 
map backgrounds by changing the background file name. The file format is an American Standard Code for 
Information Interchange (ASCII) latitude/longitude file each user can edit. The model also has the capability to 
generate output for Geographic Information System (GIS) applications. 
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3. Automated HYSPLIT Configuration 

The Tcl/Tk GUI available in the HYSPLIT source code is convenient for managing and configuring single use 
simulations. However, NWS MLB and SMG forecasters require automated trajectory and/or concentration guidance 
to be made on a regularly scheduled basis, so they requested a Linux version of the software that can be automated 
for real-time applications. A Linux version of HYSPLIT would also allow NWS MLB and SMG to make use of 
meteorological data already located on their local servers. The AMU contacted and obtained HYSPLIT version 4.8 
from the model developer, Dr. Roland Draxler of NOAA. This version is flexible enough to be run on a PC or 
compiled on a Linux system. 

The AMU compiled the source code using both the GNU Fortran compiler and Portland Group Fortran (PGF) 
compiler at the request of NWS MLB and SMG. This allows the forecasters the flexibility of using the code on 
several different Linux platforms and, more importantly, allows compatibility with other software programs built 
with each compiler. The NWS MLB version of HYSPLIT and several preprocessors were built using the GNU 
Fortran compiler; while the SMG version of the code was compiled with PGF compiler due to their local 
configuration of WRF to ARL converter as described in Section 4.1. The AMU tested and configured the newly 
compiled versions of HYSPLIT on the local AMU Linux cluster before building all model executables at NWS 
MLB and SMG. 

Once the HYSPLIT source code was configured the resulting executables were tested with the Tcl/Tk GUI to 
ensure the model would run successfully on Linux. At the request of NWS MLB, the AMU configured several 
scripts to automatically download several NCEP model products, convert the meteorological grids into ARL format, 
run the model from multiple latitudes, longitudes, and elevations, and post-process the data to create output 
graphics. The scripts were configured to run on a scheduled basis using cron jobs, which enable users to execute 
commands or scripts automatically. The scripts control the automatic downloads of both the NCEP RUC and NAM 
model data. HYSPLIT was also configured to run using the GFS guidance, but due to its coarse horizontal 
resolution, the forecasters requested the NAM and RUC to be their preferred input data. The cron jobs reference 
central parameter files for each NCEP product that contain many of the necessary trajectory and concentration 
HYSPLIT model runtime variables. This allows forecasters to quickly make changes to the model configuration 
without having to edit the automated scripts. 

The cron jobs configured for NWS MLB launch three different automated processes on a scheduled basis. The 
scheduled times are a compromise between when the NAM and RUC model guidance becomes available on the 
NCEP FTP server and computational limitations on the NWS MLB Linux machine. Mr. Peter Blottman, of NWS 
MLB, developed a master script for the NAM and RUC guidance to oversee each step in the automated HYSPLIT 
system. The master script checks for successful completion of each step contained in the cron jobs before the next. 
This ensures the next step will not run if there are significant delays in download times or missing data at any step 
during the process. The master script also ensures that two or more implementations of HYSPLIT do not interfere 
with each other by creating a temporary lock file at the start of the script. The lock file is deleted at the end of the 
master script after all the steps in the cron are completed successfully. The details of each step run via the cron jobs 
are described below. 

3.1 Gribmaster Program 

The first step in the automated HYSPLIT system is to obtain the necessary NCEP forecast products. These files 
are downloaded from the NCEP FTP server using the gribmaster utility program written by Dr. Robert Rozumalski 
at the NOAA NWS Science and Training Research Center (STRC). The gribmaster program is a Perl utility that 
allows users to automatically download multiple forecast times of operational model data including the GFS, NAM 
and RUC GRIB files from the NCEP FTP server. In addition, the gribmaster program has an option to convert the 
newer NCEP GRIB version 2 files into the older GRIB version 1 files. This option was activated since the 
HYSPLIT GRIB to ARL conversion utilities only accept GRIB version 1 files. 

The AMU configured the gribmaster utility scripts to automatically download the operational NAM 12-km and 
RUC 13-km forecast products. For shorter-range (< 12-hour) HYSPLIT simulations, NWS MLB prefers to use the 
RUC to drive HYSPLIT since it is updated hourly and assimilates many types of observations, including surface and 
aircraft data, at a higher temporal resolution than either the NAM or GFS. The RUC is used by the NWS, aviation, 
and severe weather community for short-range forecasts. The RUC was designed to provide numerical forecast 
guidance for a 12-hour period (Benjamin et al. 1998). It assimilates observations every hour at the surface and aloft 
and is run every hour, which is the highest frequency of all forecast models at NCEP. The AMU configured the 
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RUC gribmaster script to obtain products every 3 hours, with guidance every 3 hours out to 12 hours. After the 
GRIB version 2 files are downloaded from the NCEP FTP server the gribmaster utility converts them into GRIB 
version 1 format. 

For medium range, (< 48-hour) HYSPLIT simulations, NWS MLB prefers the NAM 12-km model over the 
operational GFS products due to the better spatial resolution. The NAM model produces an 84-hour forecast every 6 
hours, or four times per day. The AMU configured the NAM gribmaster utility program to download the NAM 12- 
km model guidance four times per day, each at 3 -hour forecast intervals out to 48 hours. The NAM gribmaster 
utility also converts GRIB version 2 files into the older version 1 format. 

NCEP has the NAM model output available in smaller subsets, called tiles, instead of the entire Continental 
United States (CONUS) domain (Figure 1). For the NWS MLB configuration, the AMU modified the original NAM 
gribmaster program to download NAM tile 16, which is centered over much of Florida (Figure 1). This option saves 
NWS MLB computational download time along with significant hard drive space on their Linux system. There are 
no such tiles available for the operational RUC guidance, so the entire CONUS domain is used. 


^6 


^8 

■j; 

50 

r"T — — 

If 

51 

■■/■'V 

52 

j- / 

Tyf'yK 

53 

.5^ 

37 

38 

\i ' 7 

■HQ 

m 

■ 

B 


:J : ; 

17’' k * 

^5 A:.' * 

■ / 

28 

29 

i .30 

3 1 

■' V ’ r ~‘"! '-2 

: 32 

v . ■ • . ,V/ : ■«' 

33 -rr 

i \A 

.j 

•• f* ;7 . 'a 

i "" 

H 

B 

19 

20 

}s %- 

\ 21 V 

\ j 

'j 

;■ j 

-A.23- - 

2*1: ■ 

■■■ 

% 26 

! 

27 

10 

1 1 

12, A 

V 13 1 

j 

v H 1 

- ’ 13* ...J 

7 f 

x * 


\ A 

\ * s 

I 

1 

17 

18 

01 

02 

03 


i 

05 \ 

06 " 

A^L.. i 

i oA ; 

^ 1 1 j 

P9- ,: r 


Figure 1. The NCEP operational NAM 12-km model domain. The output is available in the smaller tiles numbered 
from 1 to 54. 

3.2 GRIB to HYSPLIT ARL 

After downloading the NAM and RUC guidance from the NCEP FTP server, the cron jobs convert the GRIB 
files into ARL format. The NAM 12-km tiles are converted into the ARL format using the etal2arl Fortran program 
available in the HYSPLIT source code. The etal2arl program does not allow users to convert multiple NAM- 12km 
tiles into a single ARL file. If users want to use multiple NAM- 12km tiles to drive HYSPLIT the scripts would have 
to be configured to create an ARL file for each NAM tile. The etal2arl program converts the GRIB files on pressure 
levels to the ARL format. Only a single model time can be converted at once, so the AMU configured a script called 
naml2arlscr to run the etal2arl program to convert each downloaded GRIB file into ARL format.. Currently, the 
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script is configured to convert NAM tile 16 so the other tiles would need to be downloaded from the NCEP FTP 
server via the gribmaster program. The etal2arl program converts the GRIB files on pressure levels to the ARL 
format. The naml2arl.scr then combines the various ARL files into a single larger file to cover the desired forecast 
period. A single 12-km ARL file containing all desired forecast times is then used to drive both the NAM HYSPLIT 
trajectory and concentration products. It is important to note the NAM analyses fields do not contain diagnostic 
surface fields such as pressure, temperature, and 10-m winds, therefore no HYSPLIT guidance is started at those 
times. However, a forecast field from the previous 6-hour NAM cycle is used to fill in the gap. 

The RUC 13-km GRIB files are converted using the ruc2arl Fortan conversion program. This code was not 
available in the default HYSPLIT source code, so the AMU obtained and modified existing code from Mr. Glenn 
Rolph at NOAA. The ruc2arl program was originally written to convert the older RUC 20-km grids into ARL 
format, so the source code array dimensions were enlarged to support the higher resolution grid configuration. It is 
important to note there is an error output message displayed routinely upon running the ruc2arl program due to 
missing mean sea level pressure in the analysis hour of the RUC GRIB files. Even though this variable is missing in 
the model output, HYSPLIT does not use mean sea level pressure in any trajectory or concentration calculations, 
rather it uses surface pressure or terrain height from the meteorological input file so the message can be ignored. 
This issue was confirmed and discussed with Dr. Draxler and is a result of a QC routine within the ruc2arl code that 
analyzes the meteorological variables in the input file. Similar to the NAM 12-km HYSPLIT guidance, the RUC 13- 
km GRIB to ARL conversion is controlled by a script, called ruc2arl.scr that combines multiple times into a single 
larger file and is used to drive the requested RUC HYSPLIT simulations. 

3.3 HYSPLIT Runtime Scripts 

The final step in the Linux cron jobs is to produce automated trajectory and/or concentration guidance using the 
ARL formatted files created in the previous step. The scripts are configured for guidance from a pre-determined 
number of source locations throughout the NWS MLB County Warning Area and can be easily tuned for any 
location. The scripts read in central parameter files that allow the user to configure all the critical HYSPLIT 
variables for the automated runs. More importantly, the parameter files allow forecasters to tune the model without 
having to edit the automated scripts in real-time. A complete description of the variables contained in each 
parameter file can be found in Appendix A. Even though there are separate parameter files for each NCEP model 
product, the options in each parameter file are identical. This allows forecasters to configure HYSPLIT differently 
for each NCEP model product. 

The runtime scripts are configured to automatically provide trajectory and concentration guidance using the 
variables defined in the parameter files. For each NCEP product, the cron jobs run the trajectory model first and then 
the concentration model afterwards. Output graphics are created for each source location and are date- stamped 
according to the NCEP model cycle. After the concentration model is run for each desired source location the ARL 
files are deleted. This ensures the directory structure is clean and ready for the next HYSPLIT scheduled runtime. 

The NCEP trajectory runtime scripts are called run_naml2_traj.scr and run rucl 3_traj.scr. Each script does 
the following: 

1 . Reads in the NAM or RUC parameter file. 

2. References latest ARL file created by earlier cron steps. 

3. Creates the baseline CONTROL file for the trajectory model. The CONTROL file contains the number of 
source locations, latitude, longitude and elevation of each source location, trajectory lengths, data 
directories, and output file names. There is no need to edit the script since it passes these variables in 
automatically from the central parameter file for each NCEP model type. A complete description of the 
variables contained in each parameter file can be found in Appendix A. 

4. Creates an advanced configuration file for trajectory and dispersion guidance. For default simulations the 
advanced file is not required; however it contains variables such as model time step, number of vertical 
levels, and trajectory end point intervals. The script passes the variables in automatically from the central 
parameter file. 

5. Runs the HYSPLIT trajectory executable ( hymodelt ) for each requested forecast time. 

6. Creates postscript graphics file, converts it to Graphics Interchange Format (GIF), and creates Google Earth 
output files using the HYSPLIT trajplot program. The postscript file is deleted after the conversion to GIF 
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format. The graphics are then named by NCEP model type, model cycle time, and forecast hour. A single 
graphics file is created no matter how many source locations are requested. This enables forecasters to view 
multiple trajectories on a single image. Figure 2 shows an example of three 12-hour forecast trajectories 
initialized in real-time by the RUC 13-km guidance on a Google Earth map. 

7. Creates GIS shapefiles if the traj shpfile option is activated in each parameter file. This is done by re- 
running the traplot program and requesting an Earth Resource Surveys Incorporated (ESRI) formatted text 
file. The text file is then converted to shapefile format using the ascii2shp conversion utility available 
within HYSPLIT. Each parameter file contains options to create trajectory shapefiles consisting of points or 
lines. The shapefiles are then named by NCEP model type, model cycle time, and forecast hour and written 
to the requested output directory. 
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03/16/2009 2230 UTC*' 
03/1^/2009 21 30 UTC 


03/1 7/2009. 0230 UTC 


■v 

03/J 7/2009 0300 UTC 


3/17/2009 0300 UTC 


y . t 03/16,2009 2000 UTC 

03/16/20091830 UTC 0 3/1 7/2009 0300 UTC 
6/2009^800 UTC -.03/17/2009 0200 UTC 
03 -S/2009 1700 UTC 03/1 7/2009 01 00 UTC Source Locatjon 


03/17/2009 0300 UTC 


03/16/20 09 1 600, U/TC , . 03/.1 6/2009 23301UTC . 

C 3 - 6/200*9 1 500 UTC- 0 3 1 6/2009 2200 UTC 

03/16/2009 2100 UTC 


_ 0'3/,1 6/2009 2030 UTC 

\ 03/16/2009 1900 UTC 

E — s, 

03/16/2009 1800 UTC 

'\03/1 6/2009 1630-UTC 

03/16/200911 53'0‘UTC 03716/2009 1500 UTC 


Figure 2. Example HYSPLIT trajectories showing 12-hour forecasts from three source locations within the NWS 
MLB forecast area. The HYSPLIT trajectories begin at 1500 UTC 16 March 2009 and are initialized using RUC 13- 
km forecast data. Users can click on each trajectory end-point to get the time, latitude, longitude, and elevation 
along each trajectory. 
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The concentration runtime scripts are called run naml 2 cone. scr and run rucl 3_conc.scr. Each script does the 
following: 

1 . Reads in NAM or RUC parameter file. 

2. References latest ARL file created in earlier cron step. 

3. Configures the model to release pollutants at user-defined source locations over a defined number of 
forecast hours. The concentration model is run separately for each source and start-time; however the 
number of iterations is controlled in the parameter file. A complete description of the variables contained in 
each parameter file can be found in Appendix A. 

4. Creates the baseline CONTROL file for the concentration model. This file contains the source location 
information, emission rate and duration, sampling grid information, number of pollutants, and deposition 
information. These variables are configured in central parameter files. 

5. Creates an advanced configuration file for the concentration guidance. The file controls model variables 
including time step, model type such as particle, puff or Gaussian, number of particles or puffs released, 
time splitting information, and a chemistry model option. All variables are set in the parameter file and are 
passed into the script at runtime. 

6. Runs the HYSPLIT concentration executable ( hymodelc ) for each source location and forecast hour. 

7. Creates postscript graphics files, converts it to GIF format, and creates Google Earth output files using the 
HYSPLIT concplot program. The postscript file is deleted after the conversion to GIF format. The 
graphical output files are named by the NCEP model type, source identification, model cycle time, and 
forecast hour. At the request of NWS MLB, separate graphics files are created for each source location 
because multiple plumes could overlap and cause confusion to the forecaster. The Google Earth files allow 
for animated guidance products depending on the output interval defined in the parameter file. Figure 3 
shows output from the HYSPLIT concentration model for a single source release west of Melbourne, FL of 
a pollutant initialized in real-time by the RUC 13-km guidance on a Google Earth map. 

8. Creates GIS shapefiles if the conc shpfile option is activated in each parameter file. This is done by 
running the concplot program and requesting an ESRI-formatted text file. The text file is then converted to 
shapefile format using the asciHshp conversion software available within HYSPLIT. The parameter file 
contains options to create concentration shapefiles consisting of polygons or lines. The shapefiles are then 
named by NCEP model type, source identification, model cycle time, and forecast hour. It is important to 
note that activating the concentration shapefile option can take a long time to process depending on the 
resolution and number of sampling times requested in the parameter files. 

9. Removes the HYSPLIT ARL file for each NCEP model type. Note this step is not done in the trajectory 
runtime scripts as a new ARL file will be created at the next download time, which is controlled by the 
automated cron jobs. 
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Figure 3. Example HYSPLIT concentration output showing a potential pollutant release of time-integrated air 
concentrations averaged between 2000 and 2100 UTC March 2000. The concentrations are shown in the 0 to 50-m 
layer above ground as defined in the parameter files. HYSPLIT was initialized using real-time RUC 13-km forecast 
data at 1500 UTC 16 March 2009. The contour levels plotted in this example are 1 0" 10 , 10' 11 , and 10' 12 mass units 
and are available in the Google Earth legend. Even though this image represents a static concentration plot during a 
12 hour forecast, animated concentration output is available by using the animation toolbar in Google Earth. 
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4. HYSPLIT WRF EMS Configuration 

NWS MLB and SMG run locally-configured versions of the WRF model. The model consists of two dynamical 
cores, the Advanced Research WRF (ARW) and Nonhydrostatic Mesoscale Model (NMM). The ARW core was 
developed primarily at the National Center for Atmospheric Research, while the NMM was developed at NCEP. 
The NWS MLB and SMG forecasters run the WRF model using the EMS system (Rozumalski 2006). The EMS is a 
complete, full-physics, and numerical weather prediction package that incorporates dynamical cores from both the 
WRF ARW and NMM models into a single end-to-end forecasting system. The EMS includes all the pre-compiled 
binaries necessary to run all pre- and post-processing programs so no source code compilation is required on each 
user’s platform. More details on the WRF EMS can be found at htt p://strc.comet.ucar.edu/wrf/index.htm . 

The local versions of the WRF EMS are run at finer spatial and vertical resolution than either the RUC or NAM 
model guidance. Additionally, the NWS MLB version of WRF is initialized with the Local Analysis and Prediction 
System (LAPS; McGinley 1995). LAPS is a 3D weather analysis program that integrates data from meteorological 
observations into a gridded meteorological background field such as the RUC model. These locally tuned WRF 
models can provide input into HYSPLIT that may be a better depiction of the mesoscale environments than the 
coarser resolution NCEP guidance. 

Older versions of HYSPLIT relied on output from operational synoptic and mesoscale models provided by 
NCEP; however, recent modifications have allowed users the flexibility of using guidance from WRF simulations. 
Recent source code was made available to convert ARW model output in Network Common Data Form (NetCDF) 
format to ARL format; however no such code exists to use NMM output. 

4.1 SMG 

SMG uses the ARW within the EMS package so the converter recently made available in HYSPLIT could be 
used without modifications. However, the ARW to ARL source code required NetCDF libraries to be built on 
Linux. During the installation on the local AMU Linux cluster, the AMU noted that NetCDF could only be compiled 
using the PGF compiler so the ARW to ARL converter was also built using the PGF compiler. However, during the 
HYSPLIT installation at SMG, it was discovered that the NetCDF libraries could be built on their Linux system with 
the GNU Fortran compiler while the HYSPLIT source code could be built with the PGF compiler. At the request of 
SMG, the PGF complier was used to compile their version of HYSPLIT. Therefore, it is important to note the SMG 
version of HYSPLIT differs from the NWS MLB version of the model due to the Fortran compilers used in each 
case. If NWS MLB decides to use the ARW to ARL converter in the future it may require NetCDF and the 
converter to be built with the PGF compiler as was found on the AMU Linux cluster. 

The ARW to ARL source code directly converts model native output files in NetCDF format to ARL format. 
The code was made flexible enough to allow any horizontal and vertical grid configuration. The only limitation is 
users can only convert one ARW grid configuration and model output time upon execution. To alleviate this 
restraint the ARW ARL files can be combined together to link multiple model output times. The AMU developed a 
script, called arw2arl.scr , to convert ARW output into ARL format and to combine multiple times into one file. The 
AMU also tested the newly compiled source code to ensure the HYSPLIT model would run correctly off example 
ARW model output. HYSPLIT and the ARW converter were recently installed and configured on the SMG Linux 
system and were tested successfully. 

4.2 NWS MLB 

NWS MLB is running the NMM version of WRF within the EMS. The HYSPLIT conversion program already 
developed for the ARW could not be used to convert NWS MLB WRF output due to the different horizontal and 
vertical coordinate systems between the ARW and NMM. The AMU investigated whether existing source code 
could be modified or if a new program needed to be developed to ingest the locally configured NMM output. After 
investigation, the AMU decided to modify the existing Fortran code for converting the NAM operational output 
(nams2arl) into ARL format. The modifications included changing the horizontal and vertical grid arrays within the 
code to match the operational grid configuration used at NWS MLB. It is important to note the new nams2arl 
program converts NMM GRIB output as opposed to the native model NetCDF output. This requires the users to run 
the WRF post processing program (WRFPOST), the last step within the EMS, in order to convert the model output 
files from NetCDF into GRIB format. This is an important difference between the NWS MLB and SMG WRF to 
HYSPLIT conversion programs. 
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The NMM to ARL conversion program allows for any desired WRF horizontal grid configuration; however, the 
model sigma levels are currently hardwired to the current operational NWS MLB vertical grid structure. It would be 
beneficial to have the program read in the vertical grid configuration at runtime since it would allow forecasters to 
use different NMM vertical grid configurations without having to recompile the HYSPLIT source code. However, 
this was not done in the current task because it would require the WRFPOST program contained within the EMS to 
write out the individual sigma levels within the output GRIB files. Since the EMS contains pre-compiled binaries, 
changes to the WRF NMM output structure were not possible without going through the EMS developer. The AMU 
contacted Dr. Rozumalski about making this change in future versions of the EMS; however this was not addressed 
before the task ended. 

For compatibility with the EMS, the new nams2arl routine required several changes to the EMS Perl scripts and 
configuration files. The changes included: 

1 . Activating the WRF native level model output in the wrfcntrl.parm file. This file controls what variables 
are written to the GRIB files within the EMS. The nams2arl code needs model variables directly on the 
native vertical grid coordinate system. By default, the EMS version of WRFPOST outputs variables on 
pressure levels, so it was critical that native model levels be output as well. The nams2arl program required 
the following variables on model levels: pressure, height, temperature, turbulent kinetic energy, specific or 
relative humidity, u and v wind components, and a vertical wind component. 

2. Keeping the de-staggered Awakara E-grid WRF GRIB, or referred to as native model output files in the 
EMS scripts. By default, the wrfpost.pl EMS script deletes these files after interpolation to a normal 
latitude/longitude grid for ingestion into other meteorological software programs. This required 
deactivating the line in the wrfpost.pl script within the EMS system that deletes the de-staggered GRIB 
files. 

Once the above changes were made to the EMS, the AMU wrote a script to convert NMM model output into 
ARL format similar to one developed for SMG. The AMU also configured a HYSPLIT runtime script and parameter 
file similar to those developed to handle the operational NCEP guidance. 

Even though an NMM to HYSPLIT utility program was developed in this task, the operational implementation 
has not yet been employed due to an issue within the EMS. During initial testing of the new conversion program, the 
AMU discovered a bug with the treatment of NMM vertical levels within the EMS. The model output vertical levels 
were written in a top-down approach instead of the normal ground-up configuration contained in the operational 
version of WRFPOST run at NCEP. The files created by the WRF EMS postprocessor caused a runtime error when 
ingested into HYSPLIT. The AMU debugged this issue with Dr. Rozumalski and it has since been corrected in the 
latest version of the EMS, which was released in February 2009. The new EMS was tested and configured on the 
AMU Linux cluster; but due to other customer priorities it has not yet been installed at NWS MLB. It is important to 
note the newly created NMM to HYSPLIT conversion program developed in this task only works with the newest 
version of the WRF EMS. Based on discussions with NWS MLB, it was deemed beyond the scope of this task to 
modify the older version of the EMS to be compatible with HYSPLIT. However, NWS MLB is expected to upgrade 
their EMS version in the near future. Once the upgrade is complete it will be important that the changes outlined 
above are made. 
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5. Summary 

The first goal of this task was to obtain and install the latest version of HYSPLIT on a Linux system that ingests 
routine NCEP model products. The second goal of this task was to configure a utility program to convert WRF EMS 
output into HYSPLIT binary format. Both of these objectives were met successfully. NWS MLB uses HYSPLIT to 
provide trajectory and concentration guidance during the significant release of harmful chemicals, radiation, and 
smoke from fires and/or toxic plumes into the atmosphere. SMG forecasters have also expressed interest in using 
HYSPLIT to assist with Shuttle FR, in the event of low-level smoke plumes, cumulus clouds associated with smoke 
plumes, and high-level trajectory forecasts for thunderstorms. NWS MLB and SMG rely on the PC-based version of 
HYSPLIT through the use of a GUI; however it is difficult for forecasters to run and manage in an operational 
environment. To alleviate this difficulty the AMU configured the latest version of HYSPLIT on Linux that ingests 
several gridded meteorological forecast products. 

The AMU wrote several scripts that are scheduled to run through Linux cron capability to produce automated 
HYSPLIT trajectory and concentration guidance from the RUC and NAM models. The scripts are configured to 
automatically download the NAM and RUC guidance, convert the meteorological grids into ARL format, run the 
model from several latitude/longitude sites, and post-process the data to create output graphics. The scripts reference 
parameter files for each NCEP product that contain many of the necessary trajectory and concentration HYSPLIT 
runtime variables. This allows forecasters to make changes to the model configuration without having to edit the 
automated scripts. 

NWS MLB and SMG run local versions of the WRF model within the EMS that are run at finer spatial and 
vertical resolution than provided by the NCEP model guidance. This can provide input into HYSPLIT that may 
better predict the unique mesoscale environment in Florida and southern California to support Space Shuttle 
landings. SMG uses the ARW within the EMS package so a converter recently made available in HYSPLIT could 
be used without modifications. On the other hand, NWS MLB is running the WRF NMM within the EMS. No 
similar conversion program for HYSPLIT existed for this dynamic core. Therefore, the AMU created a conversion 
program for NMM output by modifying existing source code used to convert the NAM products. Even though this 
program was developed and tested as outlined in the task directive, the operational implementation was not 
employed at NWS MLB. This was due to an error in the way the WRF model vertical levels were written to the 
GRIB files in the EMS. The bug has since been corrected in the latest version of the EMS, which will be employed 
soon at NWS MLB. This will allow forecasters to automatically provide trajectory and concentration guidance on a 
scheduled basis using either the NCEP products or high-resolution local WRF output and, therefore, provide timely 
information on hazardous conditions to their customers. 

Forecasters at NWS MLB and SMG have expressed interest in a follow on task to develop a new GUI that 
would allow users to more easily make any needed adjustments to the central parameter files developed in this task 
through an intuitive interface. The GUI would support adjustments to the automated scheduled guidance from the 
NAM, RUC, and/or WRF products through changes to the parameter files that were developed in this task. 
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Appendix A 

Local HYSPLIT Directory Structure and Parameter Files 

This appendix provides a brief description of the local HYSPLIT directory structure and parameter files. The 
installation on Linux includes various sub-directories and files needed for model execution. It is important to note 
the directory structure as outlined below is similar to the default HYSPLIT configuration, but has additional 
directories due to the automated scripts. There are several differences between the directories and files at NWS 
MLB and SMG due to the local configuration of WRF at each location. For a more complete description of the 
directory structure and code, refer to Draxler 2007. Many of the directories also contain Readme files that provide a 
description of the contents beyond what is outlined here: 

ARW - contains the processed ARW ARL files for each scheduled FIYSPLIT run. This directory also contains the 
arw2arlscr conversion scripts. 

bdy files - files for controlling gridded land use, roughness length, and terrain data. The horizontal resolution of each 
product is currently one degree; however the source code is flexible enough to allow users to input their own data. 
The files are read in at model execution so there is no need to recompile the source code. 

browser - contains Tcl/Tk GUI interface code for the HyperText Markup Language (HTML) browser and is used to 
display help files. Note these files are referenced by the main HYSPLIT GUI used to make single use simulations. 

cluster - includes several sub-directories for configuring HYSPLIT on a computing cluster. 

cmaq- source code to convert HYSPLIT concentration output into the Community Modeling and Analysis System 
(CMAQ) model format. 

csource -source code used by the default HYSPLIT GUI in the guidcode directory. Several archived satellite 
databases are also included in the directory. 

data2arl - includes example programs to convert gridded meteorological data, including the ARW, in various 
formats to the ARL format. 

document- includes various readme and help files for configuring HYSPLIT, including a copy of Draxler 2007. 

examples- contains several example scripts and batch files that can be used to create automated simulations. 

exec - contains all the HYSPLIT executable programs, including the trajectory and concentration models, and 
graphical plotting programs. This directory also includes several programs for debugging and viewing the 
meteorological variables contained in each ARL file. Several of these programs were used extensively by the AMU 
in the testing phase of this task. They include: 

chkjtimes -examine the model times that are each ARL files. 

chk data -read and unpack the ARL formatted meteorological data. 

chk Jlle -examine the header and data records of an ARL file. 

chkrec —dumps the first 50 bytes of each meteorological data record. 

contour -creates a postscript file for a specified variable in the ARL file with contours and colors for a 
specified variable and time. 

cdit_miss - interpolates missing times within the ARL formatted meteorological file. 

graphics -source code and libraries for the trajectory and concentration postscript programs. This directory contains 
the map background files for all postscript plots, which is defined in the parameter files. A description of how to 
configure a customized background map can be found in the readme file contained in this directory. 

guidcode - contains a Tcl/Tk GUI interface for configuring and running HYSPLIT. The hysplit.tcl is the top-level 
script that will run all other program scripts. 

html -directory for main HYSPLIT GUI run from the guidcode directory. 

IMAGES- contains the graphical output images from each automated HYSPLIT simulation. This directory includes 
sub-directories to store the ARW, RUC, NAM, NMM guidance and is defined in the model parameter files. 
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library - includes various code libraries used to compile the model. 

NCEP - contains the NAM and RUC GRIB and ARL files for each scheduled HYSPLIT run. This directory also 
contains the naml2arl.scr and ruc2arl.scr conversion scripts. 

NMM - contains the processed WRF NMM files for each scheduled HYSPLIT run. This directory also contains the 
nmm2arl.scr conversion scripts. 

paramflle - contains the HYSPLIT parameter files used to edit and configure the automated HYSPLIT simulations. 
There is a separate parameter file for the RUC, NAM, and WRF guidance. The configuration of each file is exactly 
the same, but is used by each automated script to run the various model products. 

qwikcode - source code used to compile a version of HYSPLIT to make quick trajectory and concentration guidance 
using default model options. 

source - includes all the Fortran HYSPLIT source code. This directory includes the new ruc2arl.f and nams2arl.f 
conversion programs, which are used to convert the GRIB files into ARL format for the RUC and WRF NMM 
models, respectively. 

utility - contains the HYSPLIT build and configuration files, which are referenced upon compilation. 

In addition to the sub-directories outlined above, the main HYSPLIT directory is configured as the working 
directory for the trajectory and concentration runtime scripts and is controlled by the cron jobs. The directory is the 
default location for all written temporary files such as the CONTROL and advanced configuration files. The 
configuration files are created based upon the variables defined in the parameter files and are stored in the 
paramfiles directory. Each parameter file includes variables to define the trajectory, concentration, and graphical 
output options. A description of the trajectory and concentration variables contained within each parameter file is 
outlined in Table 1 and Table 2. The variable names were chosen to be similar to those outlined Draxler 2007. 

In most scenarios users should only have to edit the parameter files without having to make changes to the 
runtime scripts. However, in some contingency situations users may have to refer back to the process outlined in 
Sections 3 and 4 and make sure the needed meteorological input data is available to drive whatever type of 
simulation is needed. This is because the cron jobs and various scripts were designed to make automated HYSPLIT 
simulations at scheduled times. 


20 



Table 1 . A description of the HYSPLIT trajectory model variables contained in the RUC, NAM, and WRF 
parameter files. The parameter files are used by the automated scripts to write out the CONTROL and advanced 
configuration files needed to run the trajectory model. 


1 - - - 1 

Variable Name 

Example Value 

Description 

& General Parameters 


Necessary working directories 

workingdir 

/data/hysplit4 

Directory where HYSPLIT temporary 
files will be written 

gribdir or netcdf_dir 

/data/hysplit4/NCEP/NAM 

Directory where GRIB formatted 1 files 
for the RUC, NAM, and WRF NMM 
reside. For the ARW HYSPLIT 
simulations this variable contains the 
path to the model NetCDF files. 

datadir 

/data/hysplit4/NCEP/NAM/ 

Directory where processed RUC, NAM, 
and WRF ARL files reside. These files 
contain the meteorological grids used to 
drive the HSYPLIT simulation. It is 
important to keep the syntax of the 
default filename the same, including the 
last forward slash, as shown in the 
example as the HYSPLIT CONTROL 
file expects this syntax. The directory 
also contains the GRIB to ARL 
conversion scripts referenced by the 
cron jobs. 

imagedir 

/data/hysplit4/NCEP/I MAGES 

Directory where final output graphics 
are written, including Google Earth and 
any other requested graphics after 
successful HYSPLIT completion. 

shpdir 

/data/hysplit4/NCEP/SHPFILES 

Directory where output GIS shapefiles, 
if requested, are written. 

mapfile 

/data/hysplit4/graphics/floridamap 

ASCII text map background file for 
postscript graphics. This file can be 
edited or changed if the user desires. 

& Basic Trajectory Parameters 


Basic trajectory parameters. These 
variables are needed for each simulation 
and used to create the CONTROL file. 

traj_year 

00 

Two-digit start year of trajectory 
simulation. Only change this variable 
for contingency runs as when set to zero 
HYSPLIT defaults to the start of the 
meteorological guidance. The scripts 
will automatically process through the 
chosen times based upon the current 
meteorological input type and the 
variable num traj fhour. 

trajmonth 

00 

Start month of trajectory simulation. 
When this parameter is set to zero it 
defaults to the start month of the 
meteorological guidance. 
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Table 1 . A description of the HYSPLIT trajectory model variables contained in the RUC, NAM, and WRF 
parameter files. The parameter files are used by the automated scripts to write out the CONTROL and advanced 
configuration files needed to run the trajectory model. 

Variable Name 

Example Value 

Description 

trajday 

00 

Start day of trajectory simulation. When 
set to zero the model defaults to the 
start day of the meteorological 
guidance. 

trajhour 

00 

Start hour of trajectory simulation. 
When this parameter is set to zero the 
model uses the start hour of the 
meteorological guidance 

numtrajfhour 

4 

Number of start times for trajectory 
simulations. The automated runtime 
scripts will use this parameter to 
determine the number of trajectory 
iterations. 

fhourtraj 

03 

Start-times (hours) of trajectory 
simulation in relation to meteorological 
data files. The number of start-times is 
based upon the num traj Jhour 
parameter. 

traj source number 

10 

Number of source locations. 

trajlat 

28.02 

Latitude (degrees) of each source 
location. 

trajlon 

-80.10 

Longitude (degrees) of each source 
location. 

trajelev 

10 

Altitude (meters) of each source 
location. 

trajlength 

12 

Length of trajectory simulation (hours). 
Backward calculations are configured 
by setting the time to a negative value. 

traj_vert_motion 

0 

Vertical motion option for trajectory 
model, 0 - use vertical motion directly 
from meteorological input data, 

1 - derive from isobaric data, 

2 - derive from isentropic data, 

3 - constant density , 4 - constant 
internal sigma coordinate, 5 - compute 
from velocity divergence , 6 - correct 
the vertical velocity when mapped from 
quasi-horizontal ETA surfaces. 

traj mode l top 

10000 

Sets the vertical limit (meters) of 
HYSPLIT model calculations. 

traj _num_grids 

1 

Number of simultaneous input 
meteorological grids. For automated 
runs this option will not change. 

trajout 

mlb.traj 

Output trajectory file. The filename is 
used internally within the scripts, but is 
not used to name output graphics. 
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Table 1. A description of the HYSPLIT trajectory model variables contained in the RUC, NAM, and WRF 
parameter files. The parameter files are used by the automated scripts to write out the CONTROL and advanced 
configuration files needed to run the trajectory model. 
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Table 1. A description of the HYSPLIT trajectory model variables contained in the RUC, NAM, and WRF 
parameter files. The parameter files are used by the automated scripts to write out the CONTROL and advanced 
configuration files needed to run the trajectory model. 
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Table 2. A description of the HYSPLIT concentration model variables contained in the RUC, NAM, and WRF 
parameter files. The parameter files are used by the automated scripts to write out the CONTROL and advanced 
configuration files needed to run the concentration model. 


Variable Name 


Example Value 


Description 


& General Parameters 


workingdir 


gribdir or netcdfdir 


/data/hysplit4 


/data/hysplit4/NCEP/NAM 


data dir 


/data/hysplit4/NCEP/NAM/ 


imagedir 


shpdir 


map_file 


& Basic Concentration Parameters 


/data/hysplit4/NCEP/IMAGES 


/data/hysplit4/NCEP/SHPFILES 


/data/hysplit4/graphics/floridamap 


conc_year 


cone month 


Necessary working directories 


Directory where all temporary files will 
be written. 


Directory where GRIB formatted 1 
files for the RUC, NAM, and WRF 
NMM reside. For the ARW HYSPLIT 
simulations this variable contains the 
ath to the model NetCDF files. 


Directory where processed RUC, 

NAM, and WRF ARL files reside. 

These files contain the meteorological 
grids used to drive the HSYPLIT 
simulation. It is important to keep the 
syntax of the default filename the same, 
including the last forward slash, as 
shown in the example as the HYSPLIT 
CONTROL file expects this syntax. The 
directory also contains the GRIB to 
ARL conversion scripts referenced by 

the cron jobs. 

Directory where final output graphics 
are written. These include the Google 
Earth and any other requested graphics 
after successful HYSPLIT completion. 
Directory where, output GIS shapefiles, 

if requested, are written. 

ASCII text map background file for 
postscript graphics. This file can be 

edited or changed. 

Basic trajectory parameters. These 
variables are needed for each 
simulation and used to create the model 

CONTROL file. 

Two-digit start year of concentration 
simulation. When this parameter is set 
to zero HYSPLIT defaults to the start 
year of the meteorological guidance. 

The scripts will automatically process 
through the chosen times based upon 
the current meteorological guidance 

and the variable num cone fhour . 

Start day of concentration simulation. 
When this parameter is set to zero the 
model defaults to the start month of the 
meteorological guidance. 
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Table 2. A description of the HYSPLIT concentration model variables contained in the RUC, NAM, and WRF 
parameter files. The parameter files are used by the automated scripts to write out the CONTROL and advanced 
configuration files needed to run the concentration model. 


Variable Name 


Example Value 


Description 


cone day 


cone hour 


num cone fhour 


fhour cone 


concsourcenumber 

cone lat 

cone Ion 

concelev 
cone id 



cone vert motion 


conc_model_top 


concnumgrids 


cone out 


cdump 


Start day of concentration simulation. 
When this parameter is set to zero the 
model defaults to the start day of the 

meteorological guidance. 

Start hour of concentration simulation. 
When this parameter is set to zero the 
model defaults to the start hour of the 

meteorological guidance. 

Number of start times for concentration 
simulations. The automated runtime 
scripts will use this parameter to 
determine the number of concentration 
interations. 


Start-times (hours) of concentration 
simulation in relation to meteorological 
data files. The number of start-times is 

based upon num cone Jhour. 

Number of source locations. This value 
controls how many output images are 
created per time, since the source 
locations are plotted independently. 

Latitude of each source location 

Longitude of each source location. 
Altitude (meters) of each source 

location. 

Source name used to identify output 
graphics. 


Length of concentration (hours). 
Backward calculations are configured 
by setting the time to a negative value. 


Vertical motion option for 
concentration model, 0 - use vertical 
motion from meteorological input data, 

1 - derive from isobaric data, 

2 - derive from isentropic data, , 

3 - constant density , 4 - constant 
internal sigma coordinate, 5 - compute 
from velocity divergence , 6 - correct 
the vertical velocity when mapped from 

quasi-horizontal ETA surfaces. 

Sets the vertical limit (meters) of 

HYSPLIT model calculations. 

Number of simultaneous input 
meteorological grids. For automated 

runs this option will not change. 

Output concentration file. The filename 
is used internally within the scripts, but 
is not used to name output graphics. 
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Table 2. A description of the HYSPLIT concentration model variables contained in the RUC, NAM, and WRF 
parameter files. The parameter files are used by the automated scripts to write out the CONTROL and advanced 
configuration flies needed to run the concentration model. 
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Table 2. A description of the HYSPLIT concentration model variables contained in the RUC, NAM, and WRF 
parameter files. The parameter files are used by the automated scripts to write out the CONTROL and advanced 
configuration files needed to run the concentration model. 


Variable Name 


Example Value 


Description 


num vert lev 


vertlevheight 


samp_start_year 


The number of vertical levels in 
concentration sampling grid including 
the surface. 


Height of each level (meters) AGL for 
concentration output. It is important to 
keep the number of levels equal to 
numvertlev. Multiple values can be 

input here, 

Two-digit sampling start year. Each 
grid may have a different start-time. 
Only use this setting for contingency 


samp_start_month 


samp start day 


samp_start_hour 


sampstartmin 


sampend_year 


sampendmonth 


sampendday 


sampendhour 


sampendmin 


samplingtype 


sampling hour 

sampling min 

& Deposition Parameters 


Sampling start month. Each grid may 
have a different start-time. Only use 

this setting for contingency runs. 

Sampling start day. Each grid may have 
a different start- time. Only use this 

setting for contingency runs. 

Sampling start hour. Each grid may 
have a different start-time. Only use 
this setting for contingency runs. 


Sampling start minute. Each grid may 
have a different start-time. Only use 

this setting for contingency runs. 

Two-digit sampling end year. Each grid 
may have a different start-time. Set to 

zero to default to total run length. 

Sampling end month. Each grid may 
have a different start-time. Set to zero 

to default to total run length. 

Sampling end day. Each grid may have 
a different start-time. Set to zero to 
default to total run length. 


Sampling end hour. Each grid may 
have a different start-time. Set to zero 

to default to total run length. 

Sampling end minute. Each grid may 
have a different start-time. Set to zero 

to default to total run length. 

Type of concentration sampling used in 
calculations, 0 - averaging, 

1 - snapshot or instantaneous value at 
each interval, 2 - maximum 
concentration observed during 
sampling hour and/or sampling min. 


Sampling hour interval. 

Sampling minute interval. 

Deposition parameters for all emitted 
ollutants. 
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Table 2. A description of the HYSPLIT concentration model variables contained in the RUC, NAM, and WRF 
parameter files. The parameter flies are used by the automated scripts to write out the CONTROL and advanced 
configuration files needed to run the concentration model. 


Variable Name 


Example Value 


Description 


num_pollutants_dep 


part_diameter 

partdensity 

partshape 

dep_velocity 


pol_mol_weight 

0.0 

sfc reac 

0.0 

diff ratio 

0.0 

eff henry con 

0.0 

act_henry_con 

0.0 

in_cloud_rem 

0.0 

bel_c!oud_rem 

0.0 

raddecay 

0.0 

pollresp 

0.0 


& Advanced Concentration 
Parameters 


cone tratio 


Number of pollutants depositing. The 
number of depositing pollutants must 
be equal to the number of pollutants 
released in the simulation. 


Particle diameter (pm). A value of zero 

will treat the pollutant as a gas. 

Particle density (g/cc). A value of zero 

will treat the pollutant as a gas. 

Particle shape. A value of zero will 

treat the pollutant as a gas. 

Dry deposition velocity (m/s). The 
velocity can also be calculated by the 
model using the resistance method 
which requires setting the next four 
variables. Refer to Draxler 2007 for 

more details. 

Pollutant molecular weight 
(Gram/Mole' 


Surface reactivity ratio. 

Diffusivity ratio. 


Effective Henry’s Constant. 

Actual Henry’s constant for wet 

removal. 

In-cloud removal (L/L). This parameter 
is defined as a ratio of the pollutant in 

air to rain measured at the ground. 

Below-cloud removal coefficient (1/s). 
This parameter is defined through a 

removal time constant. 

Radioactive decay half-life (days). A 
non-zero value initiates decay process 
of both airborne and deposited 

pollutants. 

Pollutant re-suspension factor (1/m). A 
non-zero value causes deposited 
pollutants to be re-emitted based upon 
soil conditions, wind, and particle type. 
Advanced concentration parameters. 
These variables are not necessary for 
each simulation; since default options 
will be used. These variables are used 
to create the CONC.CFG advanced file 

at runtime. 

Defines fraction of grid cell that puff or 
particles are allowed to transit in one 
advection time step. Reducing this 
value will reduce the time step and 
increase computation time. 
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Table 2. A description of the HYSPLIT concentration model variables contained in the RUC, NAM, and WRF 
parameter files. The parameter files are used by the automated scripts to write out the CONTROL and advanced 
configuration files needed to run the concentration model. 
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Table 2. A description of the HYSPLIT concentration model variables contained in the RUC, NAM, and WRF 
parameter files. The parameter files are used by the automated scripts to write out the CONTROL and advanced 
configuration files needed to run the concentration model. 



31 






Table 2. A description of the HYSPLIT concentration model variables contained in the RUC, NAM, and WRF 
parameter files. The parameter files are used by the automated scripts to write out the CONTROL and advanced 
configuration flies needed to run the concentration model. 


Variable Name 


Example Value 


Description 


cone dxf 


conc_dyf 


cone dzf 


& Concentration Graphic Options 


cone zoom 


cone arcview 


bottom 


J * top 


contours 


display 


conc_map_proj 


cone mult 


cone circle 


conc_shp_shpfile 
conc shp type 


4+3+2+ 1 


polygons 


West-east grid factor for offsetting 
meteorological grid in ensemble 

simulation. 

North-south grid factor for offsetting 
meteorological grid in ensemble 

simulation. 

Vertical grid factor for offsetting 
meteorological grid in ensemble 
simulation. 


Graphics options for output 
concentration plots. 


Zoom level (%) for postscript graphics. 
Arcview output, 0 - none, 1 - log 10, 

2 - use explicit values, 3 - Google 

Earth. 

Bottom view layer (meters) of each 
lot. 


Top view layer (meters) of plot. 

Contour type of plot, 

0 - dynamic/exponential, 

1 - fixed/exponential, 

2 - dynamic/linear, 3 - fixed/linear, 

4 - fixed values set by user. 

Values for contour option. The 
example here shows 4 contour values. 
Option to display concentrations by, 

1- individual vertical levels, or 

2 - average levels together. 

Map projection for graphics file, 

0- Auto, 1 - Polar, 2 - Lambert 

Conformal, 3 - Mercator. 

Concentration multiplier. 

Number of 80 km rings to draw around 
each source location. This option is not 

used for Google Earth output. 

Option to create GIS shapefile format, 

0 - no, 1 - yes. 

Type of GIS shapefile plot to create, 
olygons or lines. 
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List of Acronyms 


3D 

Three-Dimensional 

AGL 

Above Ground Level 

AMU 

Applied Meteorology Unit 

ARL 

Advanced Research Laboratory 

ASCII American Standard Code for 

Information Interchange 

CMAQ 

System 

Community Modeling and Analysis 

EMS 

Environmental Modeling System 

ERSI 

Earth Resource Surveys Incorporated 

FR 

Flight Rules 

FTP 

File Transfer Protocol 

GFS 

Global Forecast System 

GIS 

Geographic Information System 

GRIB 

Gridded Binary 

GUI 

Graphical User Interface 

HTML 

HyperText Markup Language 

HYSPLIT Hybrid Single-Particle Lagrangian 

Integrated Trajectory 

LAPS 

Local Analysis and Prediction System 

MLB 

Melbourne, FL 


MSL 

Mean Sea Level 

NAM 

North American Mesoscale 

NCEP 

Prediction 

National Centers for Environmental 

NetCDF 

Network Common Data Form 

NOAA 

Administration 

National Oceanic and Atmospheric 

NWS 

National Weather Service 

QC 

Quality Control 

PC 

Personal Computer 

PGF 

Portland Group Fortran 

POR 

Period of Record 

RUC 

Rapid Update Cycle 

SMG 

Spaceflight Meteorology Group 

STRC 

Science Training and Research Center 

TCL 

Tool Command Language 

UTC 

WRF 

Universal Coordinated Time 
Weather Research and Forecasting 

WRFPOST 

WRF Post-processing 
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